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Abstract. We derive an effective action for contacts between superconducting terminals 
with arbitrary transmission distribution of conducting channels. In the case of normal- 
superconducting (NS) contacts we evaluate interaction correction to Andreev conductance and 
demonstrate a close relation between Coulomb effects and shot noise in these systems. In the 
case of superconducting (SS) contacts we derive the electron-electron interaction correction to 
the Josephson current. At T = both corrections are found to vanish for fully transparent NS 
and SS contacts indicating the absence of Coulomb effects in this limit. 



1. Introduction 

Low energy electron transport across the interface between normal metals and superconductors 
(NS) is provided by the mechanism of Andreev reflection pQ. This mechanism involves conversion 
of a subgap quasiparticle entering the superconductor from the normal metal into a Cooper pair 
together with simultaneous creation of a hole that goes back into the normal metal. Each such 
act of electron-hole reflection corresponds to transferring twice the electron charge e* = 2e across 
the NS interface and results, e.g., in non-zero conductance of the system at subgap energies [2j. 

Andreev reflection is also responsible for dc Josephson effect in superconducting weak links 
without tunnel barriers. Suffering Andreev reflections at both NS interfaces, quasiparticles with 
energies below the superconducting gap are effectively "trapped" inside the junction forming a 
discrete set of levels which can be tuned by passing the supercurrent across the system [3]. At the 
same time, these subgap Andreev levels themselves contribute to the supercurrent [SIIIIEIE] thus 
making the behavior of superconducting point contacts and SNS junctions in many respects 
different from that of tunnel barriers. 

Note that the above theories remain applicable if one can neglect Coulomb effects. In small- 
size superconducting contacts, however, such effects can be important and should in general be 
taken into account. A lot is known about interplay between fluctuations and charging effects 
in superconducting tunnel barriers [7]. Here we examine the properties of superconducting 
junctions going beyond the tunneling limit. Below we will demonstrate that Coulomb blockade 
in such junctions weakens with increasing barrier transmissions and eventually disappears in 
the limit of fully open superconducting contacts. We will also argue that in superconducting 
systems - similarly to normal contacts [8] - there also exists a direct relation between Coulomb 
effects and current fluctuations. 
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Figure 1. The system under consideration. Two big metallic reservoirs - one superconducting 
(S) and another one either normal or superconducting (M=N,S) - are connected by a short 
normal conductor. 

2. The model and effective action 

As it is shown in Fig. 1, we will consider big metallic reservoirs one of which is superconducting 
while another one could be either normal or superconducting. These two reservoirs are connected 
by a rather short normal bridge (conductor) with arbitrary transmission distribution T n of its 
conducting modes and normal state conductance Gjy = 1/-Rjv = (e 2 / 71 ") J2n^n- Both phase and 
energy relaxation of electrons may occur only in the reservoirs and not inside the conductor 
which length is assumed to be shorter than dephasing and inelastic relaxation lengths. As 
usually, Coulomb interaction between electrons in the conductor area is accounted for by some 
effective capacitance C. 

In order to analyze electron transport in the presence of interactions we will make use of an 
approach based on the effective action formalism combined with the scattering matrix technique 
QUI El El]- This approach can be conveniently generalized to superconducting systems. 
In fact, the structure of the effective action remains the same also in the supeconducting case, 
one should only replace normal propagators by 2 x 2 matrix Green functions which account for 
superconductivity, as it was done, e.g., in [71 [13| H3]. 

Following the standard procedure we express the kernel J of the evolution operator on the 
Keldysh contour in terms of a path integral over the fermionic fields which can be integrated out 
after the standard Hubbard-Stratonovich decoupling of the interacting term. Then the kernel J 
takes the form 



where (pip are fluctuating phases defined on the forward and backward parts of the Keldysh 
contour and related to fluctuating voltages Vxp across the conductor as (p\ % (t) = & Vl,2- Here 
and below we set h = 1. 

The effective action consists of two terms, S[<p] = S c [<fi] + St[<f], where 



describes charging effects and the term St[V] accounts for electron transfer between normal and 
superconducting reservoirs. It reads [14j 
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where Gm and Gs are 4x4 Green-Keldysh matrices of M- and S-electrodes which product 
implies time convolution and which anticommutator is denoted by curly brackets. In Eq. ([2]) we 



also introduced "classical" and "quantum" parts of the phase, respectively ip + = {ipx + ^2) /2 
and cp- = ip\ — tf2- 

For later purposes we also express the average current and the current-current correlator via 
the effective action as 

(l\t))=ieJv v± -^-e iS M, (4) 

where (JJ)+ = (/(t)I(f ) + I(t')I(t)). 

Let us introduce the matrix Xo [<£>+] = 1 — T n /2 + (T n /4) |Ga/, G^j |^_=o- As the action St 

vanishes for <p~(t) = one has TrlnAo = 0. Making use of this property we can identically 
transform the action Q to 

S t = ~J2Ti:]n[l + X^oX'] , (6) 

n 

where X' = 1 + {T n /A) {{Gm, &s \ ~ 2 J — -^o- Now let us separately consider NS and SS 
interfaces. 



2. 1 . NS interfaces 

At temperatures and voltages well below the superconducting gap Andreev contribution to 
the action of NS system dominates. Hence, it suffices to consider the limit of low energies 
e <C A. In this limit we can define the Andreev transmissions [2 J T n = 1^/(2 — T n ) 2 and 
Andreev conductance Ga = (2e 2 /7r)X) n ^n- Let us assume that either dimensionless Andreev 
conductance qa = ^^2 n Tn is large, qa 3> 1, or temperature is sufficiently high (though still 
smaller than A). In either case one can describe quantum dynamics of the phase variable 
<p within the quasiclassical approximation [8j [9] which amounts to expanding St in powers of 
(small) "quantum" part of the phase <p—(t). Employing the above equations and expanding St 
up to terms ~ ip 2 __ we arrive at the Andreev effective action [15] 

iS t = iS R - S 1 /, (7) 

where 

t 

iS R = -|± Jdt!<p-(f)<p + V), (8) 
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x[p A cos(2<p + (t') - 2<p+{t")) + l-p A ] (9) 

and 

(3a = E "^ (1 ~ 7;) (10) 

is the Andreev Fano factor defined in a complete analogy with the normal Fano factor 
= J2 n T n (l — T n )/ J2n T n - We observe that the action St is expressed in exactly the same form 
as that for normal conductors [8j |9] derived within the the same quasiclassical approximation 



for the phase variable <p(t). In order to observe the correspondence between the action [HI [9] and 
that defined in Eqs. ©-([9]) one only needs to interchange normal and Andreev conductances as 
well as the corresponding Fano factors 

G N o G A , Pn O /3a (11) 

and to account for an extra factor 2 in front of the phase ip+ under cos in Eq. (|9j). This extra 
factor implies doubling of the charge during Andreev reflection. 



2.2. Superconducting contacts 

Turning to superconducting contacts we assume that fluctuating phases <p±(t) are sufficiently 
small and perform regular expansion of the exact effective action in powers of these phases. 
Then we obtain 

t 

iS t = --J dt'l s ( x )v-(t') + iS R - S T , (12) 
o 

where % * s the time-independent phase difference, 



eAsin X ^ T n Ayjl - T n sm 2 ( X /2) 

IsU) = 7) / = tanh^ — . (13) 

2 n ,/l-r n sin 2 ( X /2) 2T 



defines the supercurrent across the system [161 E] and 

t t 

S R = J dt' J dt"K{t' -t")<p-(t')cp + (t"), (14) 

o o 

t t 

5/ = J dt' J dt"l{t' -t")<p-(t')<p-(t") (15) 
o o 

with both kernels 7Z(t) and X(t) being real functions. The complete expressions for these 
functions turn out to be somewhat lengthy and for this reason are not presented here. Below 
we only emphasize some of the properties of lZ(t) and I(t). 

To begin with, it is straightforward to verify that in the lowest order in barrier transmissions 
T n the result (|12|) -(|15|) reduces to the standard AES action [7] for tunnel barriers in the limit 
of small phase fluctuations. Qualitatively new features emerge in higher orders in T n being 
directly related to the presence of subgap Andreev levels ±e n (x) inside the contact. Consider, 
for instance, the kernel I(t). It can be split into three contributions of different physical origin 

l(t)=l 1 (t)+l 2 (t)+l 3 (t). (16) 

The first of these terms, Ii(t), represents the subgap contribution due to discrete Andreev states. 
The Fourier transform of this term has the form 

vrA 4 ^ f Tl sin 2 X *, x 

T 2 fl - T 

+ 



4 V l2e 2 (x)cosh 2 (e„(x)/2T) 
T 2 (l-T n )sin 4 ( X /2) 



4(x) 



l + tanh 2 (e n ( X )/2T)] [6 (w - 2e»(x)) + <^ (w + 2e n ( X ))]}. (17) 



It is obvious that this contribution is not contained in the AES action at all. The second 
term ^j(t) can be interpreted as the "interference term" between subgap Andreev levels and 



quasiparticle states above the gap. In the low temperature limit T — > the Fourier transform of 
this term differs from zero only at sufficiently high frequencies > A + e n (x)- At higher 
temperatures T > e n (x), however, X2uj vanishes only for \uj\ < A — e n {x) an d remains non-zero 
otherwise. In the limit of small barrier transmissions this term scales as X 2 oc Tn 2 and, hence, is 
not contained in the AES action either. Finally, the third term Xs(t) accounts for the contribution 
of quasiparticles with energies above the gap. In the high frequency limit w > A or for A — > 
this term reduces to the standard result for a normal conductor Z&j — > (u;/2e 2 -R/v) coth(a;/2T). 

Turning now to the function 7Z(t) in Eq. (|14[) we note that its Fourier transform can be 
represented as 1Z W = TZ'^ + ilZ'^ where both TZ'^ and TZ'^ are real functions. The function IZ'^ is 
even in oj while 7Z'^ is an odd function of to, thus implying that the function 7Z(t) is real. 

The functions TZ(t) and X(t) are not independent. For instance, the Fourier transform IZ 1 ^ is 
related to I w by means of the fluctuation-dissipation relation IZ 1 ^ = 2X U tanh(o;/2T). The two 
functions TZ' U and TZ"^ are in turn linked to each other by the causality principle: the function 
1Z{t) should vanish for t < 0. 

Finally we would like to point out that with the aid of the above Gaussian effective action 
one can easily evaluate the phase-phase correlation functions for our problem. Combining Eqs. 
(H2])-(II5D with © one finds (cf., e.g. [17]) 
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Now we will employ the above results in order to describe the effect of electron-electron 
interactions on transport properties of superconducting contacts. 



3. Coulomb blockade of Andreev reflection 

We start from NS systems. In this case in the absence of interactions we set ip + = eV and 
trivially recover the standard result I = GaV. For the current fluctuations 5I(t) from Eqs. 
d7J)-([9]) and ([5]) analogously to [8] we obtain 
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G A 



(1 " A0« coth ^ + ^ 5> ± leV) coth U ± 11 1 



2T 



(19) 



This equation fully describes current noise in NS structures at energies well below the 
superconducting gap. For eV 3> T, uj Eq. (TK)j) reduces to the result [12] while in the diffusive 
regime the correlator (|19|) matches with the semiclassical result [19| . 

Let us now turn on interactions. In this case one should add the charging term ^ to the 
action and account for phase fluctuations. Proceeding along the same lines as in [H], for gA 3> 1 
or max(T, eV) > E c = e 2 /2C we get 



I = G A V - 2e/3 j4 TIm 



H — iv^ 

2 



IV 

1 + y 



(20) 



where ^(x) is the digamma function, w = gAEc/ir 2 T + iv and v = 2eV/irT. 

The last term in Eq. (|20p is the interaction correction to the I-V curve which scales with 
Andreev Fano factor Pa in exactly the same way as the shot noise. Thus, we arrive at 



an important conclusion: interaction correction to Andreev conductance of NS structures is 
proportional to the shot noise power in such structures. This fundamental relation between 
interaction effects and shot noise goes along with that established earlier for normal conductors 
[8] extending it to superconducting systems. In both cases this relation is due to discrete nature 
of the charge carriers passing through the conductor. 

Another important observation is that the interaction correction to Andreev conductance 
defined in Eq. (|2U|) has exactly the same functional form as that for normal conductors, cf. Eq. 
(25) in [8 J. Furthermore, in a special case of diffusive systems we have Gn = Ga, Pn = Pa = 1/3 
and the only difference between the interaction corrections to the I-V curve in normal and NS 
systems is the charge doubling in the latter case. As a result, the Coulomb dip on the I-V curve 
of a diffusive NS system at any given T is exactly 2 times narrower than that in the normal case. 
We believe that this narrowing effect was detected in normal wires attached to superconducting 
electrodes [20], cf. Fig. 3c in that paper. 



4. Interaction correction to supercurrent 

Let us now turn to the electron-electron interaction correction to the equilibrium Josephson 
current (|13p . Previously such correction was analyzed in the case of Josephson tunnel barriers 
in the presence of linear Ohmic dissipation [7] . The task at hand is to investigate the interaction 
correction to the supercurrent in contacts with arbitrary transmission distribution. 

In order to evaluate the interaction correction it is necessary to go beyond the Gaussian 
effective action (|12 p - (|15p and to evaluate the higher order contribution ~ tp 3 . It is easy to 
observe that the interaction correction to the supercurrent is provided by the following non- 
Gaussian terms in the effective action: 



6{iS t ) = J J J dt 1 dt 2 dt 3 Y(t 1 ,t2,t 3 )ip-(t 1 )ip + (t 2 )<p + (t 3 ) 
+ J J J dt 1 dt 2 dUZ{t 1 ,t 2 M)^+{h)v>-{t 2 )ip-{h). (21) 
The function Y(ti,t 2 ,t 3 ) can be written as 

Y( tl ,t 2 ,t 3 ) = Jj ^^YfaM^^e-****-* (22) 

where Y(uji,lo 2 ) = Y(oj 2 ,uji). The function Z{t\,t 2 ,t 3 ) can be expressed in a similar way. 

Adding the non-Gaussian terms (|2ip to the action and employing Eq. we arrive at the 
following expression for the interaction correction 

oo oo 

f doj f duj 

5Is(x) = ie J 2^Y(w,-u){ip+(p+) u) + 2ie J —Z(0,-u)){(p+<p-) u , (23) 



where the phase-phase correlators are defined in Eq. (|18|) . 

Let us consider the first term in the right-hand side of Eq. (|23p . It is easy to see that in 
the limit of low temperatures only frequencies \u\ > A + e n (x) contribute to the integral in Eq. 
()18p for (ip + (p + ) while the contribution from the frequency interval |u;| < A + e n (x) vanishes. 
Furthermore, the leading contribution from the first term in Eq. (|23|) is picked up logarithmically 
from the interval 2A <C \uj\ <C 1/RnC where 

\U)\ 



and the function Y(u, —to) tends to a frequency independent value. 



After a straightforward but tedious calculation in the interesting frequency range oj S> A 
from Eq. ([3|) one finds 



iAsin X ^ Tjy- T n )(2 - T n sin 2 ( X /2)) , 
Y(u,-u) = 2^ T\ nr ■ 2, /0 ^3/2 F ( € n{X))- (25) 

This high-frequency term involves the factor 1 — T n , i.e. it vanishes for fully open conducting 
channels. Combining Eqs. (|24p . (|25p with (|23|) . we arrive at the expression for the supercurrent 

I(x) = ls(x) + SIs(x)- (26) 
In the limit of low temperatures the interaction correction reads 

~£ - (») si - * ? d-?:iw ( 2 - r - sin2 f ) - (27) 

where 37V = 27r/(e 2 i?7v) is the dimensionless normal state conductance of the contact. This 
result is justified as long as the Coulomb correction 6Ig (x) remains much smaller than the non- 
interacting term Is(x) (USD- Typically this condition requires the dimensionless conductance to 
be large g N > ln(l/2AR N C). 

Note that Eq. (|27|) was derived only from the first term in Eq. (|23|) . The second term in this 
equation involving the function Z(0, —oj) and the correlator (99+93-) can be treated analogously. 
It turns out to be smaller than that of the first term by the logarithmic factor ~ ln(l/2Ai?jvC)- 

Let us emphasize again an important property of the result (|27p : The interaction correction 
contains the factor 1—T n and, hence, vanishes for fully open barriers. In other words, no Coulomb 
blockade of the Josephson current is expected in fully transparent superconducting contacts. 

The expression for the interaction correction (j27|) can further be specified in the case of 
diffusive contacts. In the absence of interactions the Josephson current in such contacts follows 
from f|13|) and takes the form corresponding to the zero-temperature limit of a well known 
Kulik-Omelyanchuk formula for a short diffusive wire 



7rA x 1 + sin f 

— — — cos — In 77 

2eR N 2 1 - sin f 



J S(X) = 7T~E> cos £ In tt— T-f ■ (28) 



Including interactions and averaging (|27p with the bimodal transmission distribution 

1 

T ny /T=% 



P(T n ) oc - J-— . (29) 



one finds 



■ X , . -1 X\ , 1 + sin(x/2) 

sin — h sm — in ■ — tt — 2 

2 2) l-sin( X /2) 



(30) 



Note that the result (|27|) can formally be reproduced if one substitutes T n — > T n + 5T n into 
Eq. (fT3|) . where 

5T„ = -— In ( I ) T n (l - T n ), (31) 

and then expands the result to the first order in 5T n . Interestingly, the same transmission 
renormalization (|3ip follows from the renormalization group (RG) equations [1U|, [TT1 [T2] 

dT n T n (l-T n ) ( 1 \ 



derived for normal conductors. In order to arrive at Eq. (|31 j) one should just start the RG 
flow at e = 1/RjsfC and stop it at e = 2A. Thus, the result (j27|) can be interpreted in a very 
simple manner: Coulomb interaction provides high frequency renormalization T n + 5T n (|3ip 
of the barrier transmissions which should be substituted into the classical expression for the 
supercurrent (fl"3|) . It should be stressed, however, that the last step would by no means appear 
obvious without our rigorous derivation since the Coulomb correction to the Josephson current 
originates from the term ~ </>— in the effective action which is, of course, totally absent in 
the normal case. 

5. Summary 

In this paper we derived a general expression for the effective action of superconducting contacts 
with arbitrary transmissions of conducting channels. In the case of NS systems we described the 
interplay between Coulomb blockade and Andreev reflection and demonstrated a direct relation 
between shot noise and interaction effects in these structures. The fundamental physical reason 
behind this relation lies in discrete nature of the charge carriers - electrons and Cooper pairs - 
passing through NS interfaces. Our results allow to explain recent experimental findings [20J. 

Superconducting contacts with arbitrary channel transmissions show qualitatively new 
features as compared to the case of Josephson tunnel barriers [7] . The main physical reason for 
such differences is the presence of subgap Andreev bound states inside the system. Our results for 
the interaction correction might explain a rapid change between superconducting and insulating 
behavior recently observed [21j in comparatively short metallic wires with resistances close to 
the quantum resistance unit ~ 6.5 YSl in-between two bulk superconductors. Previously it was 
already argued [22] that such a superconductor-to-insulator crossover can be due to Coulomb 
effects. Our present results provide further quantitative arguments in favor of this conclusion. 
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